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Abstract 

The probability distribution function (PDF) tails of the zonal flow 
structure formation and the PDF tails of momentum flux by incor- 
porating effect of a shear flow in ion-temperature-gradient (ITG) tur- 
bulence are computed in the present paper. The bipolar vortex soli- 
ton (modon) is assumed to be the coherent structure responsible for 
bursty and intermittent events driving the PDF tails. It is found that 
stronger zonal flows are generated in ITG turbulence than Hasegawa- 
Mima (HM) turbulence as well as further from marginal stability. This 
suggests that although ITG turbulence has a higher level of heat flux, 
it also more likely generates stronger zonal flows, leading to a self- 
regulating system. It is also shown that shear flows can significantly 
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reduce the PDF tails of Reynolds stress and structure formation. 
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I Introduction 



In order for magnetic confinement fusion energy to be a viable source of 
energy in our daily life, there are some current problems that have to be 
resolved. Among them is the control of the anomalous transport in Toka- 
maks, which is mostly attributed to a variety of micro-instabilities such as 
the ion-temperature-gradient (ITG) mode, trapped-electron (TE) mode and 
electron-temperature-gradient (ETG) mode pQ. Micro-instabilities can how- 
ever generate secondary modes e.g. the zonal flow that in turn may reduce 
the turbulent fluctuations in the plasma [2]- [3]. Note that zonal flows are 
poloidally and toroidally symmetric flow {kg = and k\\ = 0) structures with 
a strong inhomogeneity in the radial direction {kr 7^ 0) [4J- [6J. 

On the other hand coherent structures such as blobs and streamers can 
mediate significant heat transport through the formation of rare avalanche 
like events of large amphtude, as indicated by recent numerical studies [7]- 
[9j. Such events cause the deviation of the probability distribution functions 
(PDFs) from a Gaussian profile on which the traditional mean field theory 
(such as transport coefficients) is based. In particular, PDF tails due to 
rare events of large amplitude are often found to be substantially different 
from Gaussian although PDF centers tend to be Gaussian [10]. These non- 
Gaussian PDF tails are manifestations of intermittency, caused by bursts and 
coherent structures. The characterization of these PDF tails thus requires 
a non-perturbative method. In previous works the PDF tails of momentum 
flux {R) and heat flux have been addressed [H]- [13]. It was shown that 
the PDF tail exhibits an ubiquitous exponential (~ e~^^^^^) and that the 
coefficient ^ contains all the model dependent information, giving a difference 
between the ITG and the Hasegawa-Mima (HM) models [HJ; PDF tails in 
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both cases significantly deviate from tlie Gaussian distributions for typical 
values of C,- Moreover, it was found that reversed modon speed [H] may have 
an important influence on the PDF tail [T5] . 

In this paper we predict the PDF tails of the zonal flow structure forma- 
tion and the PDF tails of momentum flux by incorporating the effect of a 
shear flow. The model used in the present study is an advanced fluid model 
for the ITG mode [16] that permits the existence of a two dimensional bipolar 
vortex soliton (modon) [13], [T7]- [T^]. Note that the ITG turbulence model 
has been successful in reproducing both experimental [20] and non-linear 
gyro- kinetic results [21j. The rare event of large amplitude represented by 
the modon can drive a zonal flow through the generalized Reynolds stress. 
The generation mechanism is studied in detail and the properties of the PDF 
tail of the zonal flow formation is studied. 

The theoretical technique used here is the so-called instanton method, a 
non-perturbative way of calculating the PDF tails. The PDF tail is first for- 
mally expressed in terms of a path integral by utilizing the Gaussian statistics 
of the forcing. An optimum path will then be associated with the creation 
of a modon (among all possible paths) and the action is evaluated using the 
saddle-point method on the effective action. The instanton is localized in 
time [22], existing during the formation of the modon. Thus, the bursty 
event can be associated with the creation of a modon. 

It is found that the PDF tail of zonal flow {4>zf) formation has a different 
exponential form (~ e~^^'^zF^ compared with the momentum flux (R) and 
heat flux (~ e~^^^^^^) and that the constant (j = 1,2) may also differ 
significantly. Furthermore the PDF tail of zonal flow formation in ITG tur- 
bulence is much stronger than that in HM turbulence. This suggests that 
stronger zonal flows are generated in ITG turbulence than in HM turbulence 
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as well as further from marginal stability. Namely, while ITG turbulence is 
a state with high level of heat flux, it also generates strong zonal flows that 
inhibit transport. This suggests that turbulence and zonal flows constitute 
a self-regulating system whereby zonal flows generated by turbulence damp 
turbulence. It is shown that the shear flow can significantly reduce the PDF 
tails of Reynolds stress and the zonal flow formation. 

The paper is organized as follows. In Sec. II the model of the ITG 
mode turbulence is presented together with preliminaries for the path-integral 
formulation for the PDF tails of structure formation. In Sec III the instanton 
solutions are calculated and the PDF tails of momentum flux and structure 
formation are estimated in Sec IV. We provide numerical results in Sec. V 
and a discussion of the results and conclusion in Sec. VI. 

II Non-perturbative calculation of structure 
formation PDF 

The ITG mode turbulence is modeled using the continuity and temperature 
equation for the ions and considering the electrons to be Boltzmann dis- 
tributed piij. The effects of parallel ion motion, magnetic shear, trapped 
particles and finite beta on the ITG modes are neglected since in previous 
works the effect of parallel ion motion on the ITG mode was shown to be 
rather weak [23j. The effect of a shear flow is incorporated in the time evolu- 
tion equations for the background fluctuations in the form of sheared velocity 
Vq. The continuity and temperature equations then become. 



d 



dy L J L J 




-[0,T,] + ^[0,n]. (2) 



Eqs. (l)-(2) are closed by using the quasi-neutrality condition. Here [A, B] = 
[dA/dx)[dB /dy) — {OA/ dy){dB / dx) is the Poisson bracket; / is a forcing; 
n = {Ln/ps)Sn/no, = {Ln/ Ps)e6(f)/Te, Ti = {Ln/ Ps)STi/Tio are the nor- 
mahzed ion particle density, the electrostatic potential and the ion tem- 
perature, respectively. In equations (1) and (2), r = Ti/T^, Ps = Cs/Qci 



where Cg = JTe/nii, Qci = eB/rriiC and z/ is coUisionality. We also define 



major radius and = r (1 + rji). The perpendicular length scale and time 
are normalized by ps and Ln/cg, respectively. The geometrical quantities are 
calculated in the strong ballooning limit {6 = 0, gi{6 = 0,R,) = 1/R where 
Qi (6) is defined by (6) = uo^enQi (0) and R is the plasma cross-sectional 
elongation) [25]- [26j, with tu^ = kyV^ = psCgky/Ln- It should be noted that 
the time evolution of the zonal flow potential {4>zf) is governed by an aver- 
aged Eq. (1) [see last line in Eq. (6)]. This is equivalent to using the electron 
density response = (f) — (0) [21] . 

We calculate the PDF tails of momentum flux and zonal flow formation 
by using the instanton method. To this end, the PDF tail is expressed in 
terms of a path integral by utilizing the Gaussian statistics of the forcing 
/ [22]. The probability distribution function of Reynolds stress Zi = R and 
zonal flow formation Z2 = (pzF (using the notation Zj for j={l,2}) can be 
deflned as 




Lf = - (dlnf/dr) (/ = {n,Tj}), r]i = Ln/Lr^, e„ = 2L„/i? where R is the 



P{Z) 



{S{Z, - Z)) 




(3) 



6 



where 

= (exp(-zA,Z,)). (4) 

The integrand can then be rewritten in the form of a path-integral as 

h. = j V(t>V4>V(t>ZFV4>zFe''^^i. (5) 

Here, the parameter j refers to the two specific cases included in the present 
study; namely j — 1 gives the PDF tail of momentum flux while j — 2 
the PDF tail of the structure formation. The angular brackets denote the 

average over the statistics of the forcing /. By using the ansatz Tj = x0 [see 
Eq. (10)], the effective action Sx. in Eq. (5) can be expressed as, 

^A, = -^/c^^a;cit0(^-(^-a,^)Vi0 + Vo(l-Vi)0 
+ (l-e„^,/3)|^-/3[0,Vi0]j 

+ (fx(fx'4){x)K{x - x')4>{x') 

+ iX2 J dt(t)zF{t)5{t) 

- ijdt4>ZF{t){^^^-r{v,Vy)). (6) 

Note that the PDF tails of momentum flux and structure formation can be 
found by calculating the value of 5"^^ at the saddle-point in the two cases 
Ai ^ oo, A2 = or Ai = 0, A2 — > oo, respectively. Here the term {v^Vy) is the 
Reynolds stress averaged over the forcing (/) and space, that generates the 
zonal flow. The flrst case Ai 00, A2 = gives the PDF tail of momentum 
flux in ITG turbulence including the interaction of a shear flow (Vq) while 
the second limit Ai = 0, A2 — > 00 gives the PDF tail of zonal flow formation. 
In Eq. (6) we have used, 

13 = 1 + T + Tx, (7) 



X = 



77,-1(1-^ + ^0) 

U-Vo + IrenQi 



(8) 



In Eq. (8), U is the modon speed [see Eq. (10)]. To obtain Eq. (6) 
we have assumed the statistics of the forcing / to be Gaussian with a short 
correlation time modeled by the delta function as 



and (/) = 0. The delta correlation in time was chosen for the simplicity of 
the analysis. In the case of a finite correlation time the non-local integral 
equations in time are needed. We will also make use of the completeness 
of the Bessel function expansion and write n{x — x') = K,o{Jo{kx) Jo{kx') + 



Ji{kx)Ji {kx') (cos 9 cos ^'+sin 9 sin 9')+J2{kx)J2{kx') (cos 29 cos 2^'+sin 29 sin 29') + 



III Instanton (saddle-point) solutions 

We have now reformulated the problem of calculating the PDF to a path- 
integral in Eq. (6). Although the path integral cannot in general be cal- 
culated exactly, an approximate value can be found in the limit Xj — > oo 
by using a saddle point method. Since a direct application of the saddle- 
point equations results in very complicated partial differential equations for 
(f) and (f), we assume that the instanton saddle-point solution is a temporally 
localized modon. Note that in Ref. [14| the properties of the modons are 
explained in detail as well as some basic statistical models for an ideal gas 
of modons with the corresponding PDFs. That is, we assume that a short 
lived non-linear vortex soliton solution exists to the system of Eqs (l)-(2), by 
assuming that the electric potential and ion temperature Tj can be written 



{fix,t)fix\t')) = 6it-t')Kix-x'), 



(9) 



...). 
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as 

(j){x,y,t) = i!{x,y -Ut)F{t), and Ti = x(j)- (10) 

Here, 

ipix^y — Ut) = ciJi(kr)(cos6 + esinO) + —rr cos9 for r < a, (11) 

ip{x^y — Ut) = C2Ki{pr){cos6 + e{r) sin6) ioi r > a. (12) 

Here Ji and Ki are the first Bessel function and tlie second modified Bessel 
function, respectively; r = y^x'^ + y'^, tanO = y'/x, y' = y — Ut, a = {Ai — 
fcMs), = (1 - e^gi -U + Vo)/P, A2 = {U + ai- Vq)/^. By matcliing tlie 
inner and outer solution at r = a we find tlie conditions Ci = —aa/ Ji{ka), 
C2 = -Ua/Ki{pa), J[{ka) / Ji{ka) = (1 + k^/p'^)/ka - kK[{pa)/pKi{pa); 
U is the velocity of the modon, and a is the size of the core region. The 
function e(r) is chosen such that the matching conditions are similar to those 
in previous previous studies [TT]- [I2]. It is important to note that when 
e = the Reynolds stress vanishes. The action Sx^ is then to be expressed 
only as an integral in time by using the conjugate variables. 



Fo = J (fx(^{x,t)Jo{kr), (13) 

Fis = J cfx(j){x,y,t)Ji{kr) sin 6, (14) 

Fic = J (fx(j){x,y,t)Ji{kr) cos9, (15) 

F2s = J d'^x(j){x,y,t)J2{kr) sin 26, (16) 

F2c = [ (fx^{x,y,t)J2{kr)cos29. (17) 



Note that the contribution from the outer solution (r > a) to Sx^ is neglected 
compared to that from the inner solution (r < a) for simplicity. This can 
be justified since the outer solution decays fast and inherently gives a minor 
contribution to the PDF tail. The action Sx^ consists of four different parts; 
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the ITG model, the forcing, the Reynolds-stress parts or the structure for- 
mation and the time evolution of the zonal flow respectively. We note that 
the spatial structure of the zonal flow is not given which would only influence 
the constant 0zf(O) in Eq. (36), acting as a normalization in Eq. (64). The 
full action including the forcing and Reynolds stress or structure formation 
terms can then be expressed in terms of F, F and the conjugate variables F, 

^A, = / dtYiiF{Fic + eFi,) + F(72F2s + 673^^0 + e74F2c) 

+ F^{^^F2s + £76^0 + e^ihc) + 78(Ac + eFu))] 

+ \k, I dt{F^ + 2{Fl + Fl) + 2{Fl + Fl)) 

+ iXiRoj dtF\t)6{t) 

+ iX2 j dt(j)zFS{t) 

- i J dt^zFi^ + RoF') (18) 
Here the coefficients are. 



7i = 




(19) 


72 = 


k 1 
= -73 = "274, 


(20) 


75 = 


~2 = -76, 


(21) 


77 = 


(3ka, 


(22) 


78 = 


vk\ 


(23) 




1 - engi(3 -U- k\U + a,) + ^0(1 + A;') 


(24) 


/ 1*1 


.} dx dy 


(25) 


f'h: — 


Ro 
Tra^ 


(26) 



Note that the effect of the shear flow {Vq) is incorporated in the parameter 
a and cti. 
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The equation of motion for the instanton is found by the variation of the 
action with respect to F, Fq, Fic, Fis, F2c and F2s, 



5F 



6S 



A2 



5F 



= -i[-7i(Fic + eFis) + (72-F2S + 673^0 + e-f4F2c) 

+ 2F(75F2s + 676^^0 + e77i^2c) + ^s{Fic + e^^i^)] 

- 2iRo(f)zFF - 2iXiRoFd{t) = 0, (27) 

= -i [-71 {Fic + eFis) + (72 As + £73-^ + 

+ 2F(75F2s + 676^^0 + e^7F2c) + -fsiFu + eFi.)] 

- 2iRo^zFF^0, (28) 
= -ie{-fsF + 'jeF^) + KoFo = 0, (29) 

= -i(7iF + 78F) + 2fi;oFie = 0, (30) 

= -ie(F + 78F) + 2«o A. = 0, (31) 



SSx- 
SFo 

SSx- 
SFu 

s s 

= -ie(74F + 77F2) + 2/^0^20 = 0, (32) 

0F2C 
S S 

^ = -^^J^F + J,F') + 2KoF2s = 0, (33) 

s s 

-j^ = <j)zF + RoF^^O, (34) 

OCPZF 

SSx^ 

S4>ZF 

5Sx2 



4>zF = 0, (35) 

. , ^ZF + A20ZFo5(t) = 0. (36) 

dcpZF 

Here the initial condition for the zonal flow is used as 0zfo = 0zf(O). The 

equation of motion for F is derived for t < using Eqs. (27)- (36) as, 

1 dP"^ 

-7?(l + e')-TTr = 771F + 3772^2 + 2773F=^-4«o0zFi^, (37) 

2 dr 

where 

771 = 72+2eM + e'7l + 78+e'78, (38) 
?72 = 7275 + 2e^7376 + e^7477, (39) 
773 = 7| + 2e^62 + e^72. (40) 
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Ordering the terms on the RHS of Eq. (37) in powers of A2 gives that the 
dominating terms are 2rjsF^ and 4Ko(j)zFF in the structure formation hmit 
(A2 —>■ 00 and Ai = 0). The contribution from the dissipation [u) to the term 
involving the time derivative of F (F) cancels out and the equation of motion 
is exactly solvable. In the hmits of A2 — > 00 and Ai = 0, Eqs (27)-(36) give 



us, 



F 
A 

C 



FVAF^ + c, 
2K0A2 



+ 

Separation of variables then leads to, 
1 



(41) 
(42) 

(43) 



In 



VF^ + Co-VC'oJFi + Q 



/Co 



(44) 



Here Co = C/A. The general solution to Eq. (44) can be written in terms of 
H as 



H{t) 



±Vct 



where 



H{t) = 



Hoe 



VF'' + Co - 
VF'^ + Co + VC'o' 



(45) 



y^^~+Co 



(46) 
(47) 



Fq + Co — -\/C^ 

The time dependent function F can now be determined by using Eq. (45) 
(47), 

, 2x/Co^e±^* 



Fit) 



Ho = 4A-2. 



(48) 
(49) 
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Eq. (48) is derived from Eq. (44) by solving for F(t) and substituting H(t) 
using Eq. (45). The factor that determines the locahzation of F{t) in time 
is given by the initial condition for the conjugate variable of the zonal flow. 
The initial condition is found by integrating Eq. (36) over the interval [—6, 0] 
{6 ^ 1) and observing that the conjugate variables mediating between the 
forcing and F vanish for t > 0. The initial condition for the zonal flow can 
be written, 

^zf{-S) + = 0. (50) 

For t < 0, we have (pzF = 0- We use this initial condition to compute 
the saddle-point action and then predict the scaling of Sx^ (as A2 00) to 
compute the PDF tail in the next section. The limit where Ai — 00 and 
Ai = is treated in a similar way as in Ref. [15] to compute the PDF tails of 
momentum flux. Following the steps in Eq. (37) - (50), but using a different 



^~ {i+%r ^^^^ 

we can determine the initial condition for the modon (F(0)) for the PDF tail 
of momentum flux. The initial condition for the zonal flow in this case is 
found from Eq. (38) or (p^p = 0. Here the zonal flow and the conjugate 
zonal flow potential {4>zf and (pzp) is evolving passively following the drift 
wave and the conjugate drift wave potential (0 and (p). 



IV The PDF tails 



The PDF tails are found by calculating the value of using Eq. (18) at 
the saddle-point in the two cases; the PDF tail of momentum flux by taking 
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into account the effect of a shear ffow (Ai oo, A2 = 0) and the PDF tail of 
structure formation of zonal flow (Ai = 0, A2 ^ 00). The integral in Eq (18) 
is divided in four parts Ci the ITG integral; C2 the forcing integral; C3 the 
momentum flux integral; C4 the zonal flow structure integral and finally C5 
represents the zonal flow evolution integral. The resulting integrals in the 
first case are; 

^A, ^ -IhXl (52) 
h = C1 + C2 + C3 + C4 + C5, (53) 

+ 24(i.U2.S?)^(i(^-i^)). (57) 

C3 = i?oi^'(0), (58) 

_ 2i?ov/C 1 

- A H,-V ^^^^ 

Ho = 4A-2, (60) 

^ - (T^- (''^ 

Note that here C4 vanishes. In determining the integrals in Eq. (18) we have 
used the locahzed modon solution in Eq. (48). By using C defined in Eq. 
(51) and A in Eq. (61), we can find the coefficient h in terms of the modon 
parameters using Eqs (19) - (25). The PDF tail of the Reynolds stress {R) 
can now be found by performing the integration over Ai in Eq. (6) using the 
saddle-point method in the same fashion as done in Ref. i.e. recall Eq. 
(6) gives P(i?) Jd\iexp{-\iR-SxA ~ / c/Ai exp{-Aii?+/iA?}. Now, the 



saddle point integral is evaluated at the maximum point Ximax = \J R/ (3/^) 
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with the result 



PiR) ~ eM-U^r% (62) 
■ ' ' (63) 



3V3^ 

This result is similar to the previous results where a similar exponential PDF 
was found. However, here the coefficient is modified by the presence of 
the shear flow {Vq). 

In the second hmit where Ai = and A2 — > oo gives the PDF tail of the 
structure formation itself, the action integral yields 

^A. ^ —hXf, (64) 

C4 = 2(l>zFo^. (65) 

Here, the parameter h is deflned above in Eq. (53) by using C4 given in Eq. 
(65) and here C3 = 0. The PDF tail can now be computed in a similar way 
as done previously with the result 

P(0zf) ~ exp{-^2{(l)ZF/(l>ZFof}: (66) 

Recall that in the computation of the PDF tails of the zonal flow structure 
formation, the zonal flow was assumed to be driven by a modon in the ITG 
turbulence since the bipolar vortex soliton (modon) was assumed to be cre- 
ated by the forcing and that F{t) = as t — >■ — oo, wc could treat Eqs 
(62) and (66) as the transition amplitudes from an initial state with no fluid 
motion to a state with with different R/ Rq or (I)zf/<I>zfo- 

Note that all physical quantities are included in the parameters and ^2, 
through the ion temperature gradient {r]i) , density gradient (e„) , temperature 
ratio (r = Ti/Te), modon size (a), modon speed (U) and wave number (k). 
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It is also important to note that — > oo (i.e. PDF vanishes) as the the 
forcing disappears (kq 0); the instanton cannot form and the PDF vanishes 
(P(i?) — > and P{4>zf) —>■ 0)- In addition when the coupling constant e — > 
the Reynolds stress vanishes {Rq — > 0) with vanishing PDF tails (-P(-R) — > 
and ^ 0). 

We note that the exact form of the exponent may depend on the temporal 
and possibly spatial correlation of the forcing (/). In the present paper, the 
forcing is chosen to be temporally delta correlated for simplicity. The coef- 
ficients may change if the spatial structure of the solution is changed, i.e. 
other non-linear solutions to Eq. (l)-(2) were used. In general, for calculat- 
ing the PDF tail a weighted sum over various coherent structures is needed. 
At present, the only known exact solution is the modon which we have as- 
sumed to be the underlying coherent structure. A generalization should be 
straightforward if other non-linear solutions were available. The case where 
the forcing is non-Gaussian (/) and other spatial coherent structures will be 
addressed in future publications. 

V Results 

We have presented a theory of the PDF tail of structure formation and how 
the PDF tail of momentum flux is modified by the presence of a shear flow. 
The exponential forms of the two PDF tails are completely different, signi- 
fying the difference in the physical interpretation. In the case of structure 
formation the PDF tails are found as ~ exp{—^2(l>'zF}^ while the momentum 
flux PDF tail ~ exp{— ^li?^/^}. In this section the parametric dependencies 
of will be studied in detail. The results will be compared with a Gaussian 
prediction and the results in the case of Hasegawa-Mima (HM) turbulence. 
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First the PDF tail of momentum flux in ITG turbulence incorporating 
the effects of shear flow is shown in Figure 1. The parameters are rji = 4.0, 
r = 0.5, e„ = 1.0, g, = 1, a = 2, U = 2.0, kq = 0.3, e = 0.1, k ^ 1.91 with 
Vo = 0.0 (blue line, solid line), Vq = 4.0 (red line, dash-dotted line), Vq = 8.0 
(black hne, dotted line), Vq = 12.0 (green line, dashed line) and Vq = 16.0 
(magenta line, thick solid line). When Vq = 0, the result recovers the previous 
finding in Ref . [15] . It is clearly shown that the PDF tail of momentum flux 
is significantly reduced if a strong shear flow is present whereas weak flow can 
increase the PDF tail. In the equations (c.f. Eq. (18)-(26)) the flow speed 
and the modon velocity comes in as a combination of the form {U — Vq), 
determining the behavior of the resulting PDF tails. When (f/— Vq) decreases, 
the PDF tail increases until it eventually decreases. This means that there 
exists a negative value {U — Vq) that gives maximum PDF tail, depending on 
all other parameters. Although due to the Galilean invariance of Eq. (l)-(2) 
we may always perform Galilean transformations on the instanton solution to 
find other instanton solutions, those solutions correspond to different ground 
states in the field theory we consider and must be discarded since they have 
nonzero velocity in the infinity past [27]. 

Second, the PDF tails of structure formation as a function of zonal flow 
potential {(pzp) are shown in Figure 2. The parameters are rji = 4.0, r = 0.5, 
e„ = 1.0, Qi = I, a = 2, Vq = 12.0, U = 2.0 {U = -5.0 dotted black 
line and dashed black line), hq = 0.3, e = 0.1 and k ^ 1.84 (ITG case), 
k = 0.81 (ITG with reversed modon speed), k ^ 1.73 (HM case) and k = 1.56 
(HM with reversed modon speed). The PDF tails of ITG turbulence (blue 
line, solid line), forced Hasegawa-Mima (HM) turbulence (red line, dash- 
dotted line), Gaussian distribution with same parameters as for the ITG 
mode turbulence (green line, thick solid line) and cases with negative modon 
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0.5 1 1.5 2 2.5 

log(R) 



Figure 1: (Color online). The momentum flux PDF tail in ITG turbulence 
with the effects of shear flow. The parameters are rji = 4.0, r = 0.5, e„ = 1.0, 
c,. = 1, a = 2, f/ = 2.0, Ko = 3 X 10^, e = 0.1, k ^ 1.91 with Vq = 0.0 (blue 
line, solid line), Vq = 4.0 (red line, dash-dotted line), Vq = 8.0 (black line, 
dotted line), Vq = 12.0 (green line, dashed line) and Vq = 16.0 (magenta line, 
thick solid line). 
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speed in ITG mode turbulence (dotted black line) and in HM turbulence 
(dashed black line) are displayed. The PDF tail in the HM model is obtained 
by setting the parameters /3 = 1 and e„ = 0.0. This is equivalent to letting 
r)i — 0.0 and r = Ti/Tg = 0.0. The PDF tails of structure formation in 
ITG turbulence is enhanced over the HM result indicating that stronger 
zonal flows are generated in ITG turbulence compared to HM turbulence. 
This has also been seen in simulations. Note however that the Gaussian 
prediction (~ exp{— ^20|f}) is larger than the PDFs of structure formation 
(~ exp{— ^20|f}) mainly due to the difference in exponential (0|j? — ^ 0|f)- 
A reversed modon speed may enhance the PDF tail. 

Third, the PDF tail of structure formation as a function of zonal flow 
potential {(f) zf) with shear flow strength (Vq) as a parameter is displayed in 
Figure 3. The resulting PDF tails are shown for Vq = 0.0 (blue line, solid 
line), Vq — 4.0 (red line, dash-dotted line), Vq — 8.0 (black line, dotted line), 
Vq — 12.0 (green hue, dashed line) and Vq — 16.0 (magenta line, thick solid 
line). The other parameters are the same as those in Figure 1. It is found 
that the PDF tail of structure formation is decreased for large shear flow 
velocity V^, whereas weak flow can increase the PDF tail. The main reason 
for this behavior is the same as explained in Figure 1. 

Fourth, the PDF tail of structure formation as a function of zonal flow 
potential {4>zf) with rji as a parameter is displayed in Figure 4. The PDF tails 
are shown for rji — 2.0 (red line, dashed line), rji — 4.0 (blue line, solid hne) 
and r]i = 6.0 (black hne, dash-dotted line). The parameters are Vq = 14.0, 
k ~ 1.91, U = 4.0 and the others are the same as in Figure 2. The PDF tail 
is significantly increased with increasing r)i, indicating that a stronger zonal 
flow is generated with increasing normalized temperature gradient rji. This 
suggests that further from marginal stability with a higher level of heat flux. 
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Figure 2: (Color online). The PDF tail of structure formation as a func- 
tion of zonal flow potential. The PDF tails ITG turbulence (blue line, solid 
line), forced Hasegawa-Mima (HM) turbulence (red line, dash-dotted line), 
Gaussian distribution with same parameters as for the ITG mode turbulence 
(green line, thick solid line) and cases with negative modon speed in ITG 
mode turbulence (dotted black line) and in HM turbulence (dashed black 
line). The parameters are rji = 4.0, r = 0.5, e„ = 1.0, gi = 1, a = 2, 
Vo = 12.0, U = 2.0 {U = -5.0 dotted black line and dashed black line), 
Ko = 3 X 10^ e = 0.1 and k ^ 1.84 (ITG case), k = 0.81 (ITG with reversed 
modon speed), k ^ 1.73 (HM case) and k = 1.56 (HM with reversed modon 
speed). 
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Figure 3: (Color online). The PDF tail of structure formation as a function 
of zonal flow potential with zonal flow strength as a paramter. The result- 
ing PDF tails are shown for Vq = 0.0 (blue hne, sohd hne), Vq = 4.0 (red 
line, dash-dotted line), Vq = 8.0 (black line, dotted hne), Vq = 12.0 (green 
line, dashed line) and Vq = 16.0 (magenta line, thick solid line). The other 
parameters are the same as in Figure 1. 
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Figure 4: (Color online). The PDF tail of structure formation as a function 
of zonal flow potential with parameter. The PDF tails are shown for 

Tji = 2.0 (red line, dashed line), r]i = 4.0 (blue line, solid line) and rji = 6.0 
(black line, dash-dotted line). The parameters are Vq = 14.0, k ^ 1.91, 
U = 4.0 and the others are as in Figure 2. 

the micro-scale ITG turbulence is more likely to generate stronger zonal 
flows, leading to a self-regulating system. That is, turbulence and zonal 
flows constitute a self-regulating system whereby zonal flows generated by 
turbulence damp turbulence 

VI Discussion and conclusions 

The PDF tail of momentum flux ~ eyip{—^i{R/ Rq)^^'^} has been obtained 
in previous studies and seems to be ubiquitous in drift wave turbulence [TT] - 
|15j . This scaling holds when the interaction with shear flow is incorporated 
as shown in this paper. The origin of this scaling is the quadratic non- 
linearity in the dynamical system (Eqs 1-2). The difference in scaling of the 
PDF tails momentum flux and structure formation comes from the change 
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in the temporal behavior of the modon (a change in the initial condition Fq). 
The spatial structure of the modon is of less importance in determining the 
PDF tails than the temporal behavior. Therefore an approximate spatial 
structure is sufficient to determine the exponential scaling whereas the time 
dependency may affect the scaling. Note that the results for momentum 
flux PDF tails are also applicable to the PDF tails of heat flux, as shown 
in Ref. [15]. The spatial structure of the flow is incorporated in the initial 
condition of the flow (pzFo- 

We note that a non-Gaussian scaling of the PDF (the exponent of R) is 
found even when the forcing is Gaussian, although the exact exponent may 
depend on the temporal and possibly spatial correlation of the forcing (/). 
In the present paper, the forcing is chosen to be temporally delta correlated 
for simplicity. The exponent may also change, if another spatial coherent 
structure is introduced i.e. another non-linear solution to the Eq. (l)-(2) is 
found. The case where the forcing is non-Gaussian (/) will be addressed in 
a future publication. 

In summary, this paper presents the first prediction of the PDF tails of 
structure formation. One of the important results is that the PDF tail of 
structure formation from ITG turbulence is significantly increased compared 
to that in HM turbulence. Zonal fiows are also shown to be more likely 
to be generated further from marginal stability, which will then regulate 
ITG turbulence, leading to a self-regulating system. Namely, while ITG 
turbulence is a state with high level of heat fiux, it also generates stronger 
zonal fiows that inhibit transport. This also suggests that stronger zonal 
flows are generated in ITG turbulence compared with ETG turbulence. It 
was also shown that shear flows can signiflcantly reduce the PDF tails of 
Reynolds stress and zonal flow formation. 
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